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Abstract 

New method of calculation of master integrals using differential equations and 
asymptotical expansion is presented. This method leads to the results exact in 
space-time dimension D having the form of the convergent power series. As an 
application of this method, we calculate the two-loop master integral for "crossed- 
triangle" topology which was previously known only up to 0(e) order. The case 
when a topology contains several master integrals is also considered. We present 
an algorithm of the term-by-term calculation of the asymptotical expansion in this 
case and analyze in detail the "crossed-box" topology with three master integrals. 



1 Introduction 

High precision of modern high-energy experiments stimulates interest to two-loop calculations. These 
calculations, in comparison with one-loop calculations, are technically much more complicated. The 
tensor reduction and integration- by-parts (IBP) method [1,2] effectively reduces the integral, corre- 
sponding to any multi-loop diagram, to the linear combination of the master integrals (Mis). The 
coefficients in this linear combination often have poles of some finite order at e = (space-time 
dimension is D = 4 — 2e) . Therefore, it is necessary to know the expansion of the master integrals 
in e, at least, a few first terms of this expansion. Successful method of calculating this expansion 
term-by-term is based on the application of the differential equations method [3,4,5,6]. The results 
for two-loop diagrams are usually expressed in terms of (generalized) harmonic polylogarithms [5,6] 
(see, however, Ref. [4]). Due to rapid increase of calculational complexity when increasing the order 
of expansion, the calculations using this method are typically done up to 0(e) terms. However, it 
may be not sufficient for the calculation of some diagrams. Therefore, it is highly desirable to cal- 
culate master integrals exactly in D. One of the most powerful methods of such calculation is the 
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application of Mellin-Barnes transformation [7,8]. However, for the case of several external invariants 
and/or many internal lines this method becomes hard to apply. 

In this paper, we present a method of obtaining the asymptotic expansion (AE) in inverse powers of 
some large scale s of multi-loop master integrals exact in space-time dimension D. Moreover, this 
expansion has a finite radius of convergence which can be determined from ODE theory, thus being 
a power series representation accessible for numerical calculations in definite kinematic region. As 
an application, we calculate series representations of two-loop diagram presented in Fig.l in two 
different kinematics. Our result can be easily expanded in e and the O(e ) term reproduces the 
known result. We also present the algorithm allowing to calculate, term-by-term, the asymptotical 
expansion of the diagrams in Fig. 3. 

Our method can be described as follows. For the calculation of the master integrals of given topology 
we differentiate them with respect to s, and express the result of differentiation (with help of inte- 
gration by parts [1,2]) in terms of Mis of the same topology and its subtopologies. Thus, we obtain, 
in general case, the linear system of first-order ODEs. We assume that the Mis of subtopologies are 
already known by means of this method or some other. When the given topology has only one MI, j, 
the system is reduced to one equation. The authors of papers [5,6] also used the differential equations 
method and searched for the solution in the form of e-expansion imposing the boundary condition at 
some finite point. Unfortunately, the results obtained by the method proposed in Refs. [5,6] quickly 
become very cumbersome with increase of the order in e. The original point of our approach is the 
choice of boundary conditions. We impose boundary condition by fixing a definite term of asymp- 
totical expansion of j, namely, the term proportional to the main asymptotics of the solution of the 
homogeneous part of the differential equation, j^. In the cases considered in this paper, choice e > 
(D = 4 — 2e) allows to fix the boundary condition by calculation of the main asymptotics of j. We cal- 
culate this asymptotics by means of standard techniques (Feynman parametrization, Mellin-Barnes 
representation, etc.). Given the asymptotical expansion of inhomogeneous term, of jh, and boundary 
condition, we find the asymptotical expansion of j. From the standard theory of ODE it follows that 
obtained power series has a finite radius of convergence. 

Our approach can be also used in the case of system of several equations. In this case, the additional 
obstacles connected with finding the homogeneous solution arise. However, we can find the AE of the 
homogeneous solution up to any finite order. The same is true for solution of initial inhomogeneous 
system. 

Thus, using our approach, one can reduce the calculation of the AE (having the form of covergent 
power series) of master integral in D dimensions to the calculation of its main asymptotics. The 
radius of convergence can be determined by the methods of ODE theory. 

The remainder of the paper is organized as follows. In Section 2, we consider the case of one MI in 
the topology. We calculate the AE of the integrals corresponding to the vertex diagrams shown in 
Fig.l. In Section 3, we describe the algorithm which allows to calculate the AE of Mis in the case 
of system of differential equations. As an example, we consider the crossed-box topology with three 
Mis depicted in Fig. 3. In the Conclusion, we discuss the main results of the present paper. 
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Fig. 1. " Crossed-triangle" topology. 
2 Topology with one master integral 



Consider inhomogeneous differential equation for the MI j with respect to external scale s 



Ts j(s) = f(s)j(s) + h(s). 

The solution of this equation is 

j(s) = j h (s) f j^(s')h{s')ds' + 3h{s)^ 1 {s )j{s ) , 



(1) 



(2) 



so 



where jh(s) = exp[/ f(s)ds\ is the solution of the corresponding homogeneous equation. 



In our approach we are interested in the asymptotical expansion for large s and assume that it is 
already known for h(s)j^ 1 (s). Then it is convenient to let Sq tend to infinity: 



J(S) =J h (s) ([ S j^(s')h(s')ds'+jH\so)j(s 0/ 



(3) 



For Mis considered in the present paper, the choice e > allows to calculate the limit Sq — > oo 
separately in the first and the second terms in Eq. (3). In this case we can integrate term- by-term 
the AE of the integrand in Eq. (3). In this stage the problem of finding of the whole AE of j is 
reduced to the calculation of its main asymptotics. In the more general case, when the limits can not 
be taken separately, the problem is reduced to the calculation of few first terms of the asymptotics 
of j. As an application of our method, we calculate the AEs of the master integral ji shown in Fig. 
1 in two different kinematics. The analytical expression corresponding to the diagram in Fig. la has 
the form ^ ^ 

Ji(s,p 2 ,q 2 ) = I r 2 l 2 (r-p) 2 (l-q) 2 (r + iy > ( 4 ) 

where s = (p + q) 2 and the usual prescription r 2 — * r 2 + iO is implied. For convenience, we assume 
that p 2 ,q 2 ,s < 0, and the final result should be understood as the analytical continuation from this 
region. For the kinematics presented in Fig. la the differential equation has the form 



dji(s,p 2 ,q 2 ) 
ds 



(3 D)y 2 2 . 1 . 2 2 
— — Ji(s,p ,q ) + —h{s,p ,q), 



(5) 
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Fig. 2. Master integrals entering Eq.(5). 

where 

h(s,p 2 , q 2 ) = -(y + p 2 )F 2 (q 2 ,p 2 , s) - (y + q 2 )F 2 (p 2 , q 2 , s) 

+ (3-^-10) (Fitfp2s) + Fi(p2 ^ s)) ^ 

+ G (p 2 > G (9 2 > • 

and y = (p-q) = (s — p 2 — q 2 )/2, A 2 = y 2 —p 2 q 2 , and Fi(pj, p 2 ,, p%) , F 2 (pl,p 2 .,pl), and G(p\) are given 
by three- and two-point integrals presented in Fig. 2 (dot corresponds to the square of denominator). 
The asymptotical expansion for F\ and F 2 can be obtained from the well-known formula for off-shell 
massless triangle with arbitrary powers of denominators [7]. This expansion and formula for G(p\) 
are presented in the Appendix. 

The homogeneous solution of Eq. (5) is 

( 2 2 \ -1/2+6 
™(- S )- 1+2e - (7) 

The choice e > allows us to take the limit separately in the first and the second terms in Eq. (3). 
Indeed, the leading asymptotics of h(s,p 2 ,q 2 ) is s e . Thus, the integrand in Eq. (3) behaves as s _1_e 
and the integral converges when its lower limit is replaced by oo. Therefore, we can confine ourselves 
to calculation of the main asymptotics of j\ that determines the limit 

(jih(so,p 2 ,q 2 )ji{s ,p 2 ,q 2 )) . (8) 

We calculate this asymptotics with the help of Feynman parametrization 

2 2 r[l + 2e] f°° Ut=i dxi 5(1 - x 12U5 )(x 13 x 24 + x 5 x 1234 )~' - i: 



' ' (Ati) d Jo (-5x3X4X5 - p 2 x 3 (x 5 xi 2 + xix 2 4) - q 2 x 4 (x 5 xi 2 + x 2 x 13 )) 1+2e 

where Xj 1; j 2i ... = x^ + Xj 2 + • • • . The region, providing the main asymptotics, is 

< x 3 ~ x 4 < xi ~ x 2 ~ x 5 ^ 1 . (10) 
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The result of integration is 



3i(s,p ,q ) 



2 2\ s— >oo 



7T 



rpc] 



(-«) 



-l+2e 



(47^6 sin (vre) sin(27re) (-g 2 ) 2e (-p 2 ) 2e ' 



(11) 



We use the following expansion 



fc + Z + 1 + 2e, fc + / + 1 + e 
fc + l + e,Z + l + e 



(p 2 /sn q y s y 
fc! Z! 



(12) 



where V [--;;;] = ^ 
obtain 



nfijnffj^ ' Substituting (12) in Eq. (3) and integrating over s, we 



3i = (-2y) 



-l+2e 



o o \ -1/2+6 



(47T) 



-D 



x T[2e]- 



7T 



;-g 2 )- 2e (-p 2 )- 2e 



esin(7re) sin(27re) 
+47r 2 cot 2 (7re)(-g 2 )" e (-p 2 )" e (-s)- 2 T[2e] £ MMl^Mr 

fc,Z=0 

r 3 oo 



fc! Z! 



k + l + 2e,k + l + l + e' 
1 + e, fc + l + e,Z + l + e 



7T° 



X 



2r[-3e]sin 3 (7re) k l mn=Q 

(~^)- 4e ' 
fc + Z + m + n + 4e 



{p 2 /s)k+ m {q 2 /s) l + n ^ 



k\ Z! ml n\ 



k + Z + 1 + 2e, fc + Z + 1 + e 
1 + e, fc + l + e,Z + l + e 



m r 



(- S )- 2 -(-g 2 )- 2 - 
fc + Z + m + n + 2e 

(_ 3 )-3 e( _ p 2 ) - e ] 

fc + Z + m + n + 3e 



m r 



m + n, m + n + 2e 
1 + e + m, l + 2e + n 

m + n, m + n — 2e 
1 + e + m, 1— 2e + n 



+ 



m + n + e,m + n — e 
— e + m, 1 + 2e + n 



| (- g )- £ (-p 2 )- £ (-g 2 )- 2e r 

fc + Z + m + n + e 



m + n — e, m + n — 3e 
— e + m, 1 — 2e + n 



i>n+m=o \ 
3er[l + e] ; 

^n+m=0 

" 3er[l + e] 



+ (p *-> q) 



(13) 



Then we make a shift fc— > fc — m, Z — > Z — n and change the summation order. The finite sums over 
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m, n can be expressed in closed forms. The final result reads 



h = (-2y) 
x |r[2e] 



/ 2 2\ -1/2+6 



—D 



y 



7T 



(-g 2 )-^(V)- 2e 



esin(7re) sin(27re) 



+ AT[2e]n 2 cot 2 (ne)(-q 2 )-^-p 2 )-^-s)^ f] MM^^r 



fc,/=0 



fc! Z! 



Z + fc + 2e,Z + fc + l + e' 
l + e,Z + l + e, fc + l + e 



+ e 2 r[- e ] 3 r[e] 5; 

fcj=0 



(p 2 /s) fc (g 2 /s)' 



fc! Z! 



x |-(- s )- 4 T 
- (-s)" 2e (-g 2 )" 2 T 



+ (-s)- 3e (-p 2 rT 



fc + Z + 4e, fc + Z + 1 + 2e ' 
1 - 3e, / + 1 + 4e, 1 + fc + e 



3-Z*2 

2e 



— k, 3e, 2e 
1 + / + 4e, 1 + 2e 



Z + A; + 1 - 2e, fc + I + 2e ' 
1 - 3e, Z + 1 - 2e, 1 + fc + e 

2e,e, fc + / + 3e, fc + Z + 1 + e 
1 - 3e, 3e, 1 + Z + 2e, 1 + k + e 



2e 3 



— fc, 3e, 2e 
-Z - fc + 2e, 1 + 2e 



+ (- s )- e (V)- e (-g 2 )- 2 T 



2e, fc + Z + 1 - e, fc + Z + e 
1 - e, 1 + Z - 2e, 1 + fc + e 



(14) 



As known, j\ is finite at D — 4. It is easy to check that e-poles in individual terms cancel and the 
finite part reproduce the AE of the well-known result of [9] . 

Let us consider the homogeneous system of the differential equations in s for all Mis appearing in 
the subtopologies of j\. The point s = oo is a regular singular point of this system. Then, from the 
theory of ODE, it follows that the above expansion has a finite radius of convergence, determined by 
the closest singularity of the coefficients of this system. In our case, it is located in the point A 2 = 0. 
Thus, the convergence region of Eq. (14) is determined by the condition 



(15) 



Now we pass to the calculation of the asymptotical expansion at large s of master integral ^(s, p 2 , q 2 ) 
corresponding to the diagram shown in Fig. lb, 



J2(s,p 2 ,g 2 ) = ji(p 2 ,s,g 2 ^ 



(16) 



The differential equation has the form 



(17) 
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where 



h(s,p 2 , q 2 ) = p 2 F 2 (q 2 ,p 2 , s) - P - (q 2 + y)F 2 (s,p 2 , q 2 ) 

+ (3D S ~ ( ^ ( 4 ) ~ 3) (P 2 + y)(FM 2 ,P 2 , s) + F^p 2 , q 2 )) (18) 
J ° :i)2 { P 2 + y)G(q 2 )G(s). 



s(D - 4) 

The homogeneous solution of Eq. (17) is 



o o\ -l/2+e 

rq z • 

y 2 



j 2h (s,p 2 , q 2 ) = (- S )- 2e (-%r 1+2e ( 1 - ) ™ • (19) 



Similar to the previous case, the choice e > allows us to take the limit separately in the first and 
the second terms in Eq. (3). The calculation of the main asymptotics of j2{s,p 2 ,q 2 ) is slightly more 
difficult than that for ji(s,p 2 ,q 2 ). Using the parametric representation we have 

2 2 _ r[l + 2e] [<*> njU dxjSjl -Exj)(x l3 x 2 A + x 5 x 12U )' 1+3e 



(47r) D 7o (-p 2 x 3 x 4 x 5 - sa;3(x5Xi2 + XxX 2A ) - q 2 Xi(x 5 x 12 + x 2 x lz )) l+2e ' 

As it is well-known, the sum in the argument of the 5-function in Eq. (20) can run over arbitrary 
subset of Xi. We choose it as J2%i — ^1235- The region, providing the main asymptotics, is 

< x 3 < xi ~ x 2 ~ Xi ~ x 5 < 1 . (21) 

After the integration over x 3 we obtain 

2 2n s^oo r[2e] (-g 2 )~ 2e f 00 rfxirfx 2 rfx 4 rfx 5 5(l - xi 25 )/3" 2e a;7 2e , ooA 
J2{s,p,q) — > 7777777 ; — / 775 , „ 7 7775 , 7 „ , (22) 



(4tt) d -s Jo ((3 + xix^iP + x 4 x 15 y- 3e 
where (5 = X\X 2 + x 2 x 5 + x 5 x±. Now we use the Mellin-Barnes representation 

ioo+0 



1 / d a r[i-3e- 3 ] r[g] 

- 3d 7 2™ (xax, r )»' 1 J 



(/3 + x 4 xi5) 1 - 3e r[l-3e] J 2ni I3 l -^~ s (x 4 x 15 ) s ' 

—ioo+0 

and take the integrals over X{. The final integration over s can be done by applying the Barnes second 
lemma (see, e.g., [10]) and results in 



2 1 (-g 2 )- 2 - r[i- e ] 3 r[2e] 

J2(S ' P ' g } ~" (4tt)" -s e (l- 2e )2r[l-3e] 3F2 



1,1,1 -e 
2 - 2e, 2 - 2e 



(24) 
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The final result for the AE of j2(s,p 2 , g 2 ), for the diagram depicted in Fig. lb, reads 



32 



(4tt) d 



\-^r 1/2+£ ( ( v)- , r[i -^ [2 i 1 3 F 2 



y 



I 



r[2e] k jz \ s , \ s 



e(l -2e) 2 r[l -3e] 
2er[l-2e] 2 ~ k ' 1 



1,1,1-e 
2 - 2e, 2 - 2e 



f ^ ( (-?T e (-<o- e <$ + (- P 2 r 2e (-<i 2 r(-sycw + (-pT 26 ^ 2 )-^-*) 2 ^? 



2er[l - 3e, 1 - 2e] 

+ (-s)- 2 *(C$ + C$) + (- P r 2e (cg + C$) + (-g 2 )- 2 ^g )] } , 



'■k,l,m,n 



n+m 



n\m\{k — m)\{l — n)\ 
1 



k + l + e 



k + I — m — n + l + 2e, k + 1 — m — n + l + e 
A; — m + 1 + e, / — n+ l + e 



(-<4fc,i,i,(A>o — ^.fc,«,o,i^>o — ^U,z,o,o) , 



(25) 
(26) 
(27) 



c k} = EE 4,i,m,n777^ r [ m + n-e,m + n- l + e, 1 - 2e - m, 1 + e - n] 



k I 



1 - e 



m=0 n=0 



7T 



'A; + / + e 

r[2c] 



sin[27re] sin[7re] T[3e] (k + l + e)k\l\ 



k + l + 3e,k + l + l + e 
k + 2e,l + l + e 



(28) 



C, 



(2) 
k,l 



k I 
m=0 n=0 

nT[2e, l-3e] 



e- 1 

n - m r m + n - 3e, m + n - 1 - e, 2e - m, 1 + e - n 

/c + / — e 



2A; 



sin[27re] (k + I - e)k\l\ 



k + l + e, k + l + l-e 
k + l-2e,l + l + e 



(29) 



C. 



(3) 
fc,Z 



k + I - 2e 



v v-, . 3me — n(l — e) — 6e 2 
: 2^ 2^^4fc,;, m ,n — - — — mT[m + n - 2e,m + n - 1 - 3e, 2e - to, 2e - nj 

m=0 n=0 

7T sin[7re] 
~ sin 2 [27re] 



r[2e,l-3e] 



-2 



(k + l - 2e)k\l\ 



kY 



k + l-e,k + l + l-2e 
k + 1 - 2e, / + 1 - 2e 



(30) 



C. 



(4) 



m=0 n=0 ^ + ^ + 2e 

+ r[l - 2e, -e, 2e]<5 m+n=0 

7T 



(e — 1) n Y[m + n, m + n — 1 + 2e, 1 — 2e — to, — e — n]5„ >0 + 



(31) 



sin[27re] 



fc + Z + 1 + e, fc + I + 1 + 2e 
/c+l + e,/ + l + e 



r[-c] 



(fc + / + 2e)fcM 



3 r 2 



1, — Z, 1 -e 
1 + fc + e, 1 + 2e 
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k I 



Ck,l — E E Ak,l,m,n 
m=0 n=0 

7T 2 r[2e] 



n(l — e) — 3me + 2e(l — e) 



k + I + 2e 



r[m + n + 2e, m + n — 1 + e, 1 — 2e — m, — 2e — n] 



sin 2 [27re] T[3e] (Jfe + Z + 2e)k\l\ 



k + I + 3e, k + I + 1 + 2e 
£; + 2e,Z + l + 2e 



(32) 



cfi = E E 



k + l 



[1 — e) mn Y[m + n — 2e,m + n — l,2e — m, —e — n]5 m+n> i 



m=0 n=0 

- mT[l - 2e, -e,2e]5 m+n=1 

7T r[l-e] -A; 
~ sin[27re] e (fc + Z)fc!Z! 



(33) 



fc + Z + 2e, A; + Z + 1 
fc + 1,Z + 1 + e 



3 r 2 



-Z,2e, 1 - e 
fc + l,l + 2e 



H+k>0 



k I 

Ck j = E E , fc '',' m 7 ' n (3me — n(l — e) — 2e(l + 2e)) m T[m + n, m + n — 1 — e, 2e — m, — 2e — n]5 m> o 



m=0 n=0 



k + l 



(34) 



£<(8) Si-k,l,m,n 

k,i - 2^ 2^ k + i 



(n(l — e) — 3me) T[m + n, m + n — 1 — e, 1 — 2e — m, 2e — n]5 m+n >i 



m=0 n=0 

+ T[l - 2e, -e, 2e] (m/2 - n e/(l - 2e))5 m+n=1 



(35) 



Again, the convergence region is determined by the condition (15). Similar to the previous case, at 
D = A we reproduce from Eq. (25) the AE of the well-known result of [9]. 



3 Topology with several master integrals 



For the case of several Mis of the same topology, when Eq.(l) turns into the system of ODE, the formal 
solution of the corresponding homogeneous system is Texp[J f(s)ds}. The asymptotical expansion of 
this solution can be usually found only term-by-term up to some finite order. Our approach can be 
also applied in this case to find any finite-order term but, usually, not the common term of the AE. 

In this Section, we generalize our approach to the case of topology with several master integrals. It 
worth noting that IBP method always leads to the differential equations for Mis with the coefficients 
being rational functions of s, and the point s = oo being the regular singular point. Thus, multiplying 
the equations by some polynomial, we can reduce them to the following vector form 

(97 

sP 1 (l/s) F + P 2 (l/s)J = H, (36) 
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Fig. 3. " Crossed-box" topology. 

where P\(x) is some polynomial, -Pi(O) = 1, and ^(x) is the polynomial matrix. The algorithm of 
finding of the AE of J term-by-term is known from the theory of ODE. Here we will present this 
algorithm in the case when the asymptotical expansion of the inhomogeneous term is power-like, 



H=£H(a) S -*, (37) 

aeB 

and the following conditions hold 

1. ojj — ctj not integer for i ^ j , 

2. ai£B, 

where CKj are the eigenvalues of ^(O). In this case, we can search for the solution of the system (36) 
in the form of power series 

J=J2C(a)s~ a , (39) 

where a is, in general, not integer, and runs over some discreet set A, restricted from below. This set, 
as well as the coefficients C(a), are to be determined. It is obvious that B C A, for convenience we 
can replace B — > A and put H(a) = if a B. Comparing the coefficient in front of equal powers of 
s, we obtain the system of recurrence relations containing C(a) , C(ct — 1) , . . . C(ct — k) with some 
finite, due to the polynomial form of Pi and P 2 , integer k 

pP 2 (a") - Pi(a")a]C(a) = H(a) , (40) 

where a~ is the shift operator, a~f(a) = f(a — 1). From Eq. (40), we can express C(a) as a linear 
combination of C(a — 1) , . . . C(a — k) and H(a), unless ct coincides with one of the eigenvalues 
of ^2(0). Therefore, the condition of series breaking from below (at sufficiently small a) results to 
A = B U {aij, ccj + 1, . . .}, where are the eigenvalues of P2(0). If a = oci, the coefficients Ck(oc) can 
not be found from Eq. (40). Due to the conditions (38), this degenerate system remains consistent, 
and we fix the coefficients Ck(a) by using the system and applying the boundary conditions. It is 
seen, that the boundary conditions are fixed by the terms of AE of J proportional to s~ ai . 

In the more general case, when there are Ins in the AE of H and/or the conditions (38) do not 
hold, we have to search the solution in the form of the series (39) with C(a) replaced by finite sum 
J2 n C (n \a)\n n s. 
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(38) 



As an example, we consider "crossed-box" topology in forward kinematics (see Fig. 3) with three 
master integrals. In particular, these Mis appear in calculations of radiative corrections to cross 
sections and of impact factors in QCD. We choose the Mis as shown in Fig. 3, i.e., 

1 f d D rd D l 

1 (2it) 2D J r 2 l 2 {r+p) 2 {l+p) 2 {l-r + q) 2 ' 1 } 
_ 2 1_ r d D rd D l 

2 ~ P (2^p J [r 2 } 2 l 2 (r + p) 2 (l + p) 2 (l-r + q) 2 ' ( ' 

1_ f d D rd D l 

3 ~ {P ' q > (2ir) 2D J rW(r + p) 2 (l + p) 2 [(l -r + q) 2 ] 2 ' ( ' 



The master integrals (41)-(43) obey the following system of differential equations 



OS 

dJ 2 



ds 



\ 2 uy 



-2yeJ x - yj 2 + [q 2 + y) J 3 + H x , 
, 3u + q 2 — p 2 



e(l - Qe)yp 



h-Uuy + 2 (p A + 4p 2 (q 2 - y) - q 2 y) e) J 2 



- e (\ 2 s - y (p 2 - q 2 ) (u - q 2 J 3 + H, 



(44) 
(45) 

(46) 



2\ 2 suy^ = yp 2 e(l - 6e) ((p 2 - 2q 2 ) (p 2 - q 2 ) - 6A 2 ) J, 

S + 2ey (2\ 2 (Ap 2 + q 2 ) - p 2 (p 2 - 2q 2 ) (p 2 - q 2 )) J 2 

+ (- (p 2 ~ <f) e P 2 y$ - 4A 2 wye - 2X 2 sue + suy(3 (p 2 + 2yj e - y)) J 3 + H 3 , 

where u = (p — q) 2 = 2p 2 + 2q 2 — s and Hi are expressed in terms of Mis of subtopologies of 
"crossed-box" topology and presented in the Appendix. In this case 



^4e2-2 



P 2 (0) 



1 -2e 
1 - 2e 



(47) 



the eigenvalues being 0:^2,3 = 4e, 1, 1 — 2e. Therefore, we have to fix the three terms of AE of 
Jii J2, J3 proportional to s~ Ae , s -1 , s~ 1+2e , for example, Ci(4e), Ci(l), C 2 (l — 2e). 

First, we calculate the asymptotics of Ji. The Feynman parametrization of Ji reads 



Ji 



nf=i tfei 5(1 - ^12345) (^12^34 + ^5^1234) 1+3e 



r[i + 2e] 

(47r) D 7o (-p 2 (x 5 a:i4X23 + x 2 x 3 x 14 + Z1Z4Z23) - 2yx 5 (x 2 x 4 - x x x 3 ) - q 2 x 5 Xi 2 x u ) 1+2e 



(48) 

Again, we assume e > 0. It is easy to see that there is no contribution oc s 4e , therefore, Ci(4e) = 0. 
The main asymptotics oc s _1 is given by the region 



^ x 5 < xi ~ 0:2 ~ a; 3 ~ x 4 ^ 1 . 



(49) 
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The integral should be treated with care. Assuming that — p 2 , — q 2 > and Im(s) ^ 0, we obtain 

, ™ Cyi) i<J7T (V)- 2 -r[e] 2 r[l-2e] 2 

Jl s ~ (4tt)^ s r[2-4e] ' 1 J 

where a = sign(Im s) = (—i/ir) ]im s ^ 00 hi(s/u). 

Now we calculate the asymptotics of J 2 . When e > the integral in Eq. (42) is infrared divergent. 
This divergence can be eliminated by subtracting from J 2 the scaleless integral 

P (2n) 2D J (r 2 ) 2 l 2 p 2 (l + p) 2 (l + q) 2 ' ( j 

This subtraction determines the analytical continuation of J 2 to the region D < 4 (e > 0). After the 
Feynman parametrization, we have 

,r[2 + 2e] 



J2 = P 2 J f[dXi6(l- X 12345 ) 



X 



Xl(x 12 x 34 + X 5 X 12U ) 



3c 



-p 2 (x 5 x 1A x 23 + x 2 x 3 x 1A + XiX A x 23 ) - 2yx 5 (x 2 x 4 - x x x 3 ) - q 2 x 5 x l2 x u ) 2+2e 
x{ 1+e (x M5 f 



(52) 



(-p 2 (x 23 x A5 + x 2 x 3 ) + 2yx 3 x 5 - q 2 x 5 x u ) 2+2e 
The main asymptotics is given by the region 

^ x 2 ~ x 3 ~ x 5 <C x\ ~ X4 ^ 1 . (53) 

The final result reads 

C 2 (l - 2e) r[2e] 7r 2 (-p 2 )~ 2e (-? 2 )" 2e 



■^2 " 



s i-2e (47r) D sin(27re) sin(7re) 



s- l+2f - - (-s)~ 1+2e . (54) 



The AE of inhomogeneous terms determines the set B, 

CO 

B — \J{k-e,k,k + e,k + 2e} . (55) 
fe=i 

As we can see, the second condition in Eq. (38) is not satisfied. However, we have checked that the 
system of recurrence relations remains consistent and, therefore, we can search for the solution in 
the form (39) with the set A defined as follows 

00 

A = \J{k~2e,k-e,k,k + e,k + 2e}. (56) 
k=i 

The set {4e, 1 + 4e, ...} is not included because Ci(4e) = 0. To save space, we do not present 
the explicit form of recurrence relations which are readily obtained from Eqs.(44)-(46). One of the 
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consequences of these relations is that 6*3(1 — 2e) = Cz(l — 2e) and, therefore, the main asymptotics 
of J 3 reads 

™ r[2e] 7r 2 (V)- 2e (-g 2 )- 2e 

^3 - 



(47r) D sin(27re) sin(7re) 



-l+2e 



-l+2e 



(57) 



4 Conclusion 



In this paper we present the approach to the calculation of the dimensionally regularized multiloop 
integrals. This approach is based on the differential equations method and asymptotical expansion. 
As a demonstration of our method, we have calculated the asymptotical expansions (14) and (25) 
of the master integral, depicted in Fig.l in two different kinematical regions. Moreover, we have 
stressed that these asymptotical expansions have a finite convergence radius determined by Eq.(15), 
thus being the power series accessible for the numerical calculations. It worth noting that this master 
integral was known only in 0(e) order so far [6] and this expansion was much more lengthier than 
our exact representation. Our results can be easily expanded up to any order of e. We have checked 
that in the limit e — > our result coincides with the well-known result of Ref. [9]. The result for 
special kinematics when one or two external momenta are on mass shell (p 2 = and/or q 2 = 0) can 
be readily obtain from Eqs. (14)-(25). The case when a topology contains several master integrals has 
also been considered. We present an algorithm of the term-by-term calculation of the asymptotical 
expansion in this case. The " crossed-box" topology with three master integrals shown in Fig. 3 is 
analyzed in detail. 



5 Appendix 



The master integrals appearing in the inhomogeneous parts of differential equations (5), (17), (44)- 
(46) are conveniently expressed via the one-loop off-shell triangle with arbitrary powers of denomi- 
nators 



J(p, v, p) 



1 



d D r 



{2ix) D J \-r 2 Y\-{p - r) 2 Y\-{q + r) 2 ]P ' 
The well-known result [7] for this integral reads 



J(p, u, p) 



-k 2 csc(-D/2 — p — p) csc(-D/2 — p — v) 



(4tt) d / 2 



T[p, v,p,D - p- v - p] 



E 

n,m=0 



(p 2 /sr(q 2 /s) 



x 



m mi 



\fj,-D/2( _ 2\D/2-fi-u 



- (-s)-^(-p 2 ) D/2_M ""r 

- (s)- u (-q 2 ) D/2 -»- p r 

+ (_ s )^/2-M-^-p r 



D — p — v — p + n + m, D/2 — p + n + m 
D/2 - p - v + 1 + n, D/2 - p - p + 1 + m 
p + n + m,D/2 — v + n + m 
D/2 - p - v + 1 + n, p + p - D/2 + 1 + m 

u + n + m,D/2 — p + n + m 
p + u- D/2 + l+n,D/2-p-p + l + m 
p + n + m, p + v + p — D/2 + n + m 
_p + v-D/2 + n + l,p + p- D/2 + 1 +m 



(58) 



(59) 
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Introducing the notation 

i T\ji + v-D/2,D/2-ti,D/2-v] 
(Ati) lj ' z Fyi, v, D — [i — v\ 

we obtain for the master integrals appearing in Eqs. (5), (17) 

F n (s } p\q 2 ) = (-l) n+1 J(l,n)J(l + n-D/2, 1,1), 
F n (p\ q 2 , s) = (-1)" +1 J(l, n)J(l, 1, 1 + n - D/2) , 

F„(g 2 ,p 2 , s) = (-1)" +1 J(l, n) J(l, 1 + n- D/2, 1) , [bl) 

The inhomogeneous part of the system (44)-(45) is expressed as follows 

H 1 = (q 2 + y)T 7 -(q 2 + 2y)T 8 , (62) 



H 2 = (g 2 -p 2 - «)y(l - 2e) 2 T 1 /4 + (g 2 - p 2 + u)y(l - 3e)(2 - 3c) (1 - 2e)T 2 /(4g 2 e) 

- y(l - 3e)(2 - 3e)(l - 2e)T 3 /(2e) + (q 2 - p 2 + u)y(l - 3e)(l - 2e)T 5 /4 
+ (P 2 - g 2 )2/(l - 3e)(l - 2e)T 6 /2 + (g 4 - p 4 + u(Ap 2 - q 2 + 3y)/2)eyT 7 

- (3(g 2 - pV + (7p 2 + q 2 )u + Auy)eyT 8 /2 , 



H 3 = (q 4 -p A - su)y(l - 2c) 2 T 1 /2 - (g 4 -p 4 + su)y(l - 3c)(2 - 3c)(l - 2e)T 2 /(g 2 e ) 
+ sy(l - 3e)(2 - 3e)(l - 2e)(sT 3 + uT 4 )/e + (p 4 - q 4 + su)y(l - 3e)(l - 2e)(T 5 + T 6 ) 
+ s v{q. — p 4 + (3p 2 — q 2 )u + 2uy)eT 7 
- y(su(9p 2 + Aq 2 + 2y) - (p 4 - q 4 )(5p 2 - Aq 2 + 2y))eT 8 , 



where 



(64) 



T 1 = G 2 (p 2 ), T 5 = F 1 (p 2 ,g 2 ,n), 

T 2 = J(l, 1) J(l, 2 - D/2)(-q 2 ) D - 3 , T 6 = F^p 2 , g 2 , s) , 

T 3 = J(l, 1) J(l, 2 - D/2)(-u 2 ) D - 3 , T 7 = F 2 (p 2 , g 2 , s) , (65) 

T 4 = J(l, 1) J(l, 2 - D/2){-s 2 ) D ^ , T 8 = F 2 (p 2 , g 2 , u) . 
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